REPORT  NO.  N ADC-90004-60 


4 


AD-A233  654 


A  MODEL  FOR  THE  THERMOPLASTIC 
ANALYSIS  OF  METAL  MATRIX  LAMINATES 


David  J.  Barrett 

Air  Vehicle  and  Crew  Systems  Technology  Department  (Code  6043) 
NAVAL  AIR  DEVELOPMENT  CENTER 
Warminster,  PA  18974-5000 


20  DECEMBER  1989 


PROGRESS  REPORT 

Period  Covering  August  1989  to  December  1989 
Task  No.  A9329-32J 
Work  Unit  No.  KG  990 
Program  Element  No.  63269F 


Approved  for  Public  Release;  Distribution  /a  Unlimited 


Prepared  for 

NAVAL  AIR  SYSTEM  COMMAND 
Department  of  the  Navy 
Washington,  DC  20361-0001 


•  i,  n 

M. 


rr ! ' 

1  i  ' 


V" 


91  4 


09  024 


NOTICES 


REPORT  NUMBERING  SYSTEM  —  The  numbering  of  technical  project  reports  issued  by  the 
Naval  Air  Development  Center  is  arranged  for  specific  identification  purposes.  Each 
number  consists  of  the  Center  acronym,  the  calendar  year  in  which  the  number  was 
assigned,  the  sequence  number  of  the  report  within  the  specific  calendar  year,  and  the 
official  2-digit  correspondence  code  of  the  Command  Officer  or  the  Functional  Department 
responsible  for  the  report.  For  example:  Report  No.  NADC-88020-60  indicates  the  twentieth 
Center  report  for  the  year  1988  and  prepared  by  the  Air  Vehicle  and  Crew  Systems 
Technology  Department.  The  numerical  codes  are  as  follows: 

CODE  OFFICE  OR  DEPARTMENT 

00  Commander,  Naval  Air  Development  Center 

01  Technical  Director,  Naval  Air  Development  Center 

05  Computer  Department 

10  Antisubmarine  Warfare  Systems  Department 

20  Tactical  Air  Systems  Department 

30  Warfare  Systems  Analysis  Department 

40  Communication  Navigation  Technology  Department 

50  Mission  Avionics  Technology  Department 

60  Air  Vehicle  &  Crew  Systems  Technology  Department 

70  Systems  &  Software  Technology  Department 

80  Engineering  Support  Group 

90  Test  &  Evaluation  Group 


PRODUCT  ENDORSEMENT  —  The  discussion  or  instructions  concerning  commercial 
products  herein  do  not  constitute  an  endorsement  by  the  Government  nor  do  they  convey 
or  imply  the  license  or  right  to  use  such  products. 


Reviewed  By: 


UNCLASSIFIED 


SECURITY  CLASS'F'CATiON  OF  TH'S  PAGE 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No  0704-0188 

la  REPORT  SECURITY  ClASSiF. CATION 

Unclassified 

lb  RESTRICTIVE  MARKINGS 

2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


Approved  for  public  release;  distribution  is  unlimited. 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NADC-90004-60 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  OF  PERFORMING  ORGANIZATION 

Air  Vehicle  and  Crew  Systems 
Technology  Department 

6b  OFFICE  SYMBOL 

(&f£licable> 

7a  NAME  OF  MONITORING  ORGANIZATION 

6c  ADDRESS  (City,  State,  and  ZIP  Code) 

NAVAL  AIR  DEVELOPMENT  CENTER 

Warminster,  PA  18974-5000 

7b  ADDRESS  (City,  State,  and  ZIP  Code) 

8a  NAME  Of-  rUNDtftG /SPCtvSC.*vi»M(j 
ORGANIZATION  NAVAL  AIR 

SYSTEMS  COMMAND 

8b  OFFICE  SYMBOL 
(if  applicable) 

AIR-931B 

9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

8c  ADDRESS  (City,  State,  and  ZIP  Code) 

Washington,  DC  20361-0001 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO  NO 

63269F  0 


WORK  UNIT 
ACCESSION  NO 

KG  990 


11  TITLE  (Include  Security  Classification) 

A  Model  For  The  Thermoplastic  Analysis  Of  Metal  Matrix  Laminates 


12  PERSONAL  AUTHOR(S) 

David  Jonn  Barrett 


13a  TYPE  OF  REPORT 

Progress 


13b  TIME  COVERED  14  DATE  OF  REPORT  (Year,  Month,  Day)  IS  PAGE  COUNT 

FROM  8/89  TO  12/89  1969  December  20 


COSATI  CODES 


SUB  GROUP 


GROUP 


04 


11 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Metal  Matrix  Composites,  Thermoplasticity,  Hot  Structures 


This  report  describes  the  development  of  a  theoretical  model  for  the  thermoplastic  analysis  of  metal  matrix 
fiber  reinforced  laminates.  The  nonlinear  analysis  is  based  on  an  incremental  form  of  classical  lamination  theory  in 
which  the  laminate  can  be  loaded  by  residual  stresses,  thermal  loads  and  edge  stress  resultants.  The  constituent 
stresses  of  the  layers  of  the  laminate  are  computed  as  the  phase  average  stresses  of  the  Composite  Cylinder 
Assemblage.  The  matrix  plasticity  is  treated  by  a  transversely  isotropic  J2  theory  with  the  exception  that  the 
response  la  found  in  an  average  sense  since  the  onset  of  plastic  flow  (von  Mises’  yield  surface)  and  the  flow  itself 
(Prager’s  Rule)  are  based  on  the  matrix  phase  average  stresses.  The  hardening  is  temperature-dependent  and 
kinematic.  The  mathematical  model  also  accounts  for  the  temperature  dependency  of  the  fiber  and  matrix  material 
properties. 


20  distribution  avA'lab  hty  of  abstract 
□  UNCLASSIFIED  UNLIMITED  SD  SAVE  AS  RPT 


22a  NAME  C  ntvful.jCT  -’.j.'.-juJ-i 

David  Jonn  Barrett 


DO  Form  1473,  JUN  86 


21  ABSTRACT  SECURITY  Cl ASS'FiCATiON 

□  DTir  ijschs  Unclassified  _ _ 


22b  TELEPHONE  (include  Area  Code)  22 c  OFFICE  Symbol 

(215)  441-3770  Code  6043 


Previous  editions  are  obsolete 

S/N  0102-LF-014-6603 


S c  CuR  ~Y  CLASSIFICATION  Of  T^S  PAGE 

UNCLASSIFIED 


i 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


NADC-90004-60 


r 


CONTENTS 


Page 


FIGURES  .  iv 

TABLES  .  v 

SYMBOLS  .  v 

INTRODUCTION .  1 

PURPOSE  .  1 

BACKGROUND  .  1 

THEORETICAL  DEVELOPMENT  .  3 

UNIDIRECTIONAL  MODEL  .  3 

Elastic  Response  .  3 

Elastic-Plastic  Response .  4 

Transversely  Isotropic  Approximation  .  5 

Elastic  Response,  Transversely  Isotropic  Approximation .  7 

Plastic  Response,  Transversely  Isotropic  Approximation .  7 

LAMINATE  MODEL .  8 

PHASE  AVERAGE  MODEL  .  11 

APPLICATIONS  .  12 

DESCRIPTION  .  12 

NUMERICAL  IMPLEMENTATION  .  12 

P130X/COPPER  .  12 

CONCLUSIONS  .  41 

REFERENCES  .  42 

Accession  For  ^ 

~NTIS  GRA&I 

DTIC  TAB  Q 

Unannounced  [“} 

i  justification _ _ 

By - 

Distribution^ _ _ 


Availability  Codes 
jAvall  and/or 
iii  Hu  st  |  Special 


N  ADC-90004-60 


FIGURES 

Figure  Page 

1  Composite  Axial  Strain  vs.  Temperature  for  Both  Temperature-Dependent 

and  Temperature-Independent  Hardening .  16 

2  Composite  Axial  Strain  vs.  Temperature .  18 

3  Composite  Transverse  Strain  vs.  Temperature  .  19 

4  Matrix  Axial  Plastic  Strain  vs.  Temperature .  20 

5  Matrix  Transverse  Plastic  Strain  vs.  Temperature .  21 

6  Matrix  Axial  Stress  vs.  Temperature  .  22 

7  Matrix  Transverse  Stress  vs.  Temperature .  23 

8  Load  Function  vs.  Temperature .  24 

9  Load  Function  vs.  Temperature .  25 

10  Composite  Axial  Strain  vs.  Temperature .  26 

11  Composite  Transverse  Strain  vs.  Temperature  .  27 

12  Matrix  Axial  Plastic  Strain  vs.  Temperature .  28 

13  Matrix  Transverse  Plastic  Strain  vs.  Temperature .  29 

14  Matrix  Axial  Stress  vs.  Temperature .  30 

15  Matrix  Transverse  Stress  vs.  Temperature .  31 

16  Composite  Axial  Strain  vs.  Temperature .  32 

17  Matrix  Axial  Plastic  Strain  vs.  Temperature .  33 

18  Matrix  Transverse  Plastic  Strain  vs.  Temperature .  34 

19  Matrix  Thickness  Plastic  Strain  vs.  Temperature  .  35 

20  Composite  Axial  Stress  vs.  Temperature .  36 

21  Composite  Transverse  Stress  vs.  Temperature .  37 

22  Matrix  Axial  Stress  vs.  Temperature .  38 

23  Matrix  Transverse  Stress  vs.  Temperature .  39 

24  Matrix  Thickness  Stress  vs.  Temperature  .  40 

iv 


NADC-90004-60 


TABLES 


Table  Page 

1  Fiber  Material  Properties  .  13 

2  Matrix  Material  Properties  .  14 

3  Numerical  Study  of  a  Pi  30x/Copper  Unidirectional  Plate  .  15 


Symbol  Definition 


SYMBOLS 


A 

B 

D 

E 

ETAN 

fvfz 

9i-  92 
h 
k 
K 


K , 


xx 


M. 


xx 


N 


xx 


KyyKxy 
’  Myy  ,  ^xy 
Nyy  ,  Nyy 


Q 

0 

s 


>ijkl 


uY,uv,  u 


x 

i° 
JX  1 


T' 

j° 

Jy, 

vf>  Km 
x,  y,  z 

aU 


z 

o 

z 


Laminate  membrane  stiffness  matrix 
Laminate  membrane-bending  stiffness  matrix 
Laminate  bending  stiffness  matrix 
Young’s  modulus 
Tangential  modulus 

Functions  representing  the  mathematical  operations  of  the  Composite  Cylinder 
Assemblage  model 

Functions  corresponding  to  the  transversely  isotropic  approximation 

Hardening  parameter 

Strength  parameter 

Matrix  of  plate  curvatures 

Normal  and  twisting  plate  curvatures 

Normal  and  twisting  plate  moments 

Normal  and  shear  stress  resultants 

Layer  stiffness  matrix  expressed  in  local  coordinates 

Layer  compliance  matrix  expressed  in  gloobal  coordinates 

Layer  compliance  matrix  expressed  in  local  coodinates 

Matrix  deviatoric  stress  tensor 

Compliance  tensor 

Plate  thickness 

Temperature 

Reference  temperature 

Plate  displacements 

Reference  surface  displacements 

Fiber  and  matrix  volume  fractions 

Global  coordinates 

Secant  thermal  expansivity 

Thermal  expansion  matrix 


NADC-90004-60 


Symbol  Definition 


r 

T 

'bey 

y°xy 

6kx'  *yy 

,°  ,  ° 
*vx’  yy 1 

e°, 

£// 

e 

©a 

°1 

°y 

°U 

4- 

v 


Equivalent  thermal  expansion  tensor 

Thermal  stress  vector  expressed  in  local  coordinates 

Thermal  stress  vector  expressed  in  global  coordinates 

Shear  plate  strain,  engineering 

Midsurface  shear  strain,  engineering 

Matrix  of  local  ply  strains 

Matrix  of  global  laminate  strains 

Normal  plate  strains 

Midsurface  normal  plate  strains 

Matrix  of  midsurface  plate  strains 

Strain  tensor 

Ply  orientation  with  respect  to  the  global  x  coordinate 

Strain  transformation  matrix 

Stress  transformation  matrix 

Stress  matrix 

Yield  stress 

Stress  tensor 

Yield  function 

Load  function 


Superscripts 

E  An  elastic  quantity 

F  A  fiber  quantity 

M  A  matrix  quantity 

P  A  plastic  quantity 

T  A  total  quantity  (elastic  plus  plastic  components) 

A  rate  quantity 

Refers  to  the  transversely  isotropic  approximation 
-1  An  inverse 


vi 


NADC-90004-60 


INTRODUCTION 


PURPOSE 

This  report  describes  the  development  of  a  thermoplasticity  model  for  the  analytical  investigation 
of  the  behavior  of  metal  matrix  fiber  reinforced  laminates.  The  nonlinear  analysis  is  based  on  an 
incremental  form  of  classical  lamination  theory  in  which  the  laminate  can  be  loaded  by  residual  stresses, 
thermal  loads,  and  edge  stress  resultants.  The  constituent  stresses  of  the  layers  of  the  laminate  are 
computed  as  the  phase  average  stresses  of  the  Composite  Cylinder  Assemblage.  The  matrix  plasticity 
is  treated  by  a  transversely  isotropic  J2  theory  with  the  exception  that  the  response  is  found  in  an 
average  sense  since  the  onset  of  plastic  flow  (von  Mises’  yield  surface)  and  the  flow  itself  (Prager’s 
Rule)  are  based  on  the  matrix  phase  average  stresses.  The  hardening  is  temperature-dependent  and 
kinematic.  The  mathematical  model  also  accounts  for  the  temperature  dependency  of  the  fiber  and 
matrix  material  properties. 

The  purpose  of  this  work  was  to  create  an  inexpensive  but  powerful  analytical  method  capable  of 
predicting  the  trends  in  the  thermomechanical  response  of  laminates  with  ductile  matrices.  Note  that 
the  method  is  not  suitable  for  the  detailed  analysis  of  structures;  for  instance,  the  model  does  not 
account  for  time-dependent  effects  or  micromechanical  damage.  But  within  the  context  of  its  limitations 
the  model  is  a  valuable  tool  for  investigating  and  interpreting  basic  materiai  behavior.  The  model  is 
therefore  highly  useful  in  the  development  of  experimental  materials  and  in  the  assessment  of  prototype 
structure. 

BACKGROUND 

Supersonic  and  transatmospheric  flight  will  cause  high  thermal  loading  of  vehicle  surfaces  and 
propulsion  systems.  Metal  matrix  composites  (MMC),  because  of  their  thermal  conductivity,  stiffness, 
strength,  and  tailorability  are  candidate  materials  of  construction  for  particularly  hot  areas  of  these 
structures.  In  the  technical  maturation  effort  of  the  National  Aerospace  Plane  program  (NASP)  the 
hardware  of  several  of  the  demonstrator  projects  has  been  constructed  of  MMC.  It  is  projected  that 
MMC  may  constitute  10%  of  the  material  of  supersonic  aircraft. 

In  contrast  to  the  desirable  properties  of  MMC  is  the  fact  that  these  materials  can  respond 
inelastically.  The  source  of  the  inelasticity  is  often  the  yielding  of  the  matrix  under  moderate  thermal 
loadings.  In  an  MMC  the  fiber  and  matrix  typically  have  greatly  different  thermal  expansivities  (the  axial 
coefficient  of  thermal  expansion  for  advanced  graphite  fibers  is  actually  negative).  Because  of  this 
mismatch  in  properties,  thermal  loadings  cause  a  buildup  of  constituent  stresses  as  the  material  tries 
to  hold  itself  together.  The  induced  stresses  are  high,  and  although  the  fiber  response  is  linear  elastic 
the  matrix  can  undergo  substantial  plastic  flow.  When  the  thermal  loads  are  removed  the  mismatch  in 
thermal  expansivity  can  cause  yielding  in  reverse  as  the  structure  cools.  The  residual  plastic  strain  that 
exists  after  a  thermal  cycle  can  grow  with  successive  thermal  exposures  leading  to  a  ratcheting  of  the 
thermal  hysteresis  loops.  The  possibility  of  high  stress  levels,  low  cycle  fatigue,  and  unacceptable  flight 
deformations  are  therefore  pertinent  design  concerns. 

A  number  of  investigations  have  been  directed  toward  studying  the  plasticity  of  MMC.  Powerful 
numerical  models  (primarily  based  on  the  Finite  Element  Method)  have  been  used  to  study  yield 
surfaces,’  nonlinear  stress-strain  response,2  and  ultimate  strength.3  Comparing  the  analytical  results 
with  experimental  data  indicates  that  plasticity  theories  can  be  used  to  predict  much  of  the  observed 
phenomena.  However,  numerical  models  by  their  detailed  nature  are  not  conducive  to  parametric 
analyses  and  can  be  expensive  to  apply  in  nonlinear  problems  that  require  iterative  solutions. 

In  order  to  create  useful  engineering  analysis  tools  a  number  of  analytical  models  have  been 
proposed.  These  models  introduce  approximations,  primarily  concerning  the  stress  fields  within  the 
constituents,  that  make  the  analysis  tractable  but  retain  the  salient  features  of  elastic-plastic  behavior. 
An  early  work  in  this  area4  considered  a  model  of  rigid  fibers  and  an  elastic-plastic  power-law  hardening 
matrix.  The  application  of  the  model  was  limited  to  the  estimate  of  transverse  stress-strain  behavior. 
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Another  more  general  approach5  6  is  known  as  the  Vanishing  Fiber  Diameter  model  (VFD)  since  it 
assumes  that  the  presence  of  the  fibers  does  not  disturb  the  shear  and  transverse  stress  fields.  By 
introducing  matrix  plasticity  this  model  was  used  to  compute  composite  yield  surfaces  and  composite 
plastic  flow.  The  VFD  performs  well  for  plane  stress  loadings  but  is  limited  in  its  ability  to  analyze  the 
effects  of  thermal  loads.  A  similar  approach7  considered  a  two-dimensional  VFD  coupled  with  a 
White-Blesseling  matrix  plasticicity  model  in  which  the  matrix  response  is  treated  by  a  number  of 
subelements  that  behave  eiastic-perfectly-plastic.  This  model  was  used  to  examine  cyclic  thermal 
loadings  but  could  not  account  for  the  ewect  of  a  full  three-dimensional  stress  state  on  the  plastic 
response  and  was  limited  to  materials  with  temperature-independent  material  properties.  Note  that  after 
an  analytical  layer  model  is  devised  it  is  a  straightforward  procedure  to  construct  an  incremental 
laminate  model  from  laminated  plate  theory  * 8 

The  work  presented  here  is  based  on  the  Phase  Average  model,10 11  which  differs  from  the  previous 
models  in  that  it  utilizes  the  Composite  Cylinder  Assemblage12  to  compute  effective  composite 
properties.  This  approach  considers  a  full  three-dimensional  stress  state  and  was  used  with  success 
in  a  previous  nonlinear  laminate  analysis13  that  examined  the  thermal  hysteresis  of  aluminum  MMC.  This 
earlier  effort  is  extended  here  to  include  temperature-dependent  properties  for  both  the  fiber  and  the 
matrix  and  a  transversely  isotropic  plasticity  model  with  temperature-dependent  hardening. 
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THEORETICAL  DEVELOPMENT 

The  formulation  starts  with  the  development  of  a  unidirectional  layer  model  for  both  elastic  and 
elastic-plastic  response  and  then  utilizes  this  model  in  the  development  of  a  general  laminate  stiffness 
relation. 

UNIDIRECTIONAL  MODEL 

ELASTIC  RESPONSE 

The  constitutive  rate  equations  for  transversely  isotropic  fibers  and  matrices  are 

1 in 


=  (2) 

where  the  S/«f/  are  fourth-order  compliance  tensors,  the  (in  are  second-order  equivalent  thermal 
expansion  tensors,  the  &,■  are  phase-average  strains,  the  on  are  phase-average  stresses,  and  T  is  the 
temperature.  The  overdcn  indicates  the  quantity  is  a  rate  terni,  the  superscripts  F  and  M  refer  to  a  fiber 
and  matrix  quantity,  and  the  superscript  E  indicates  that  the  quantity  is  elastic.  The  indices  1 ,  2,  and 
3  are  principal  material  directions  with  1  being  directed  along  the  axis  of  the  fibers.  The  equivalent 
thermal  expansion  tensors  are  computed  from 


P,f  «  x  <$+<+ aj.r  <T -T„)  (4) 

where  the  an  are  the  secant  thermal  expansivities,  TQ  is  a  reference  temperature,  and  a  comma  denotes 
differentiation  with  respect  to  the  listed  variable. 

The  effective  composite  properties  are  computed  through  the  Composite  Cylinder  Assemblage 
(CCA)  as 


(5) 

(6) 


where  the  superscript  *  indicates  a  composite  quantity,  the  functions  /,  and  /2  refer  to  the  mathematical 
operations  of  the  micromnchanical  model,  and  Vf  is  the  fiber  volume  fraction.  The  constitutive  law  for 
the  composite  is  therefore 


(7) 
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where 


Pi/  “  ckl  +  aij  +  aijJ  (T  ~T0) 


(8) 


ELASTIC-PLASTIC  RESPONSE 

It  is  assumed  here  that  the  matrix  yields  uniformly  when  the  matrix  phase  average  stress  state 
reaches  a  strain  energy  of  distortion  equal  to  the  corresponding  energy  at  yield  for  simple  tension 
Therefore,  the  onset  of  plastic  flow  is  determined  by  a  von  Mises  yield  condition  expressed  in  terms 
of  the  matrix  phase  average  stresses.  For  a  kinematic  bilinear  work  hardening  material  the  yield  function 
is  then 


4>  =  y-3k2 


(9) 


where  V  is  the  load  function  and  k  is  the  strength  parameter.  The  load  function  is  computed  from 


v  =  {Stf-h  e!fp)(S**-h  zffp) 


(10) 


where  the  S]Y  are  the  deviatoric  stresses  computed  from  the  phase  average  stresses,  h  is  a  hardening 
parameter,  and  the  superscript  P  is  used  to  designate  a  plaslic  quantity.  The  hardening  parameter  is 
computed  from 


h  =  2/3  — -  (11) 

where  E  and  £„n  are  the  matrix  elastic  and  tangential  moduli.  The  strength  parameter  is  related  to  the 
matrix  uniaxial  yield  stress  oy  through 


k2  =  219  Gy 


(12) 


Note  that  the  usual  rules  for  loading  and  unloading  apply  except  that  they  are  expressed  in  terms  of  the 
phase  average  matrix  stresses. 

The  constitutive  law  for  matrix  plasticity  is  Prager’s  Row  Rule  expressed  in  tern  is  of  the  matrix  phase 
avn  age  stresses.  The  phase  average  rr  atrix  plastic  strain  rates  are  therefore 


(13) 


,  CMP 
where  SiJkl 


is  the  fourth-order  plastic  compliance  tensor 


(14) 
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and  A 


MP 


is  the  plastic  thermal  expansion  vector 


PMP  _  h  cMrc\ir_ 
ij  -~n>T  bijkl£-kl 


MP^MP 


k’T 
h  k 


(Sfj-h  effp) 


(15) 


The  total  phase  average  strain  rates  are  found  through  the  superposition  of  the  elastic  and  plastic 
components 


+  £ 


MP 


(16) 


Substituting  equations  (2)  and  (12)  into  (16)  leads  to 


(17) 


where 


cA IT  _  cME,  cMP 
jkl  ^ijkl 


(18) 


p  iT=$i!E+$i?p  (19) 

and  the  superscript  T  refers  to  a  total  quantity. 

TRANSVERSELY  ISOTROPIC  APPROXIMATION 

In  this  work  the  CCA  is  used  to  compute  the  effective  elastic-plastic  composite  properties  from  the 
fiber  and  matrix  constitutive  laws.  However,  the  CCA  model  is  limited  to  transversely  isotropic 
constituents  whereas  the  plastic  response  can  be  fully  anisotropic.  Also  the  equivalent  thermal 
expansion  tensors  (equations  ((3),  (4),  and  (1 5)  can  also  yield  responses  more  general  than  transversely 
isotropic.  To  overcome  this  difficulty  approximations  are  introduced  to  reduce  the  response  to 
transversely  isotropic.  Let  the  approximation  be  indicated  by  the  rules  g,  and  g2  such  that 


=  *  1  (Sflj) 

(20) 

pjr=*2<p.r> 

(21) 

Pf  =  *2<P,f) 

(22) 

p.f =*2(0.f> 

(23) 

where  the  overbar  indicates  the  transversely  isotropic  approximation.  The  rules  chosen  for  g,  are 


(24a) 
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F MP  _  qMP 
*  1122  “  *  1122 

(24b) 

FMP  _  cMP 
•*  1133  ~  ■*  1122 

(24c) 

FMP  _  rMP 
•*  2222  -  ■*  2222 

(24d) 

FMP  _  qMP 
■*  2233  “  *  2233 

(24e) 

FMP  _  -MP 
■*  3333  “  •*  2222 

(240 

fMP  _  .5 

(24g) 

^ 2323  cMP  cMP 

,5  2222  ~**2233 

FMP  _  oA//> 

**  1313  “  ■*  1212 

(24  h) 

FMP  _  pMP 
■*  1212  “  **1212 

(240 

where  the  usual  symmetries  hold  and  where  all  of  the  unaddressed  terms  vanish. 

The  rules  for  gz  are 

vtr = pfr 

(25a) 

M 

xa 

(25b) 

w = 

(25c) 

with  all  other  p  terms  vanishing  (similar  results  are  obtained  for 


For  axial  loads,  thermal  loads,  or  in-plane  shear  loads  the  g,  rules  yield  exact  results.  Since  thermal 
loads  are  the  dominant  loadings  for  the  intended  application  the  approximations  are  not  compromising. 
Furthermore,  it  has  been  shown  that  for  general  loadings11  the  more  severe  approximation  of  isotropic 
plasticity  yields  tolerable  errors.  However,  as  with  all  analytical  models,  Individual  applications  must 
be  assessed  against  the  underlying  assumptions. 
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ELASTIC  RESPONSE,  TRANSVERSELY  ISOTROPIC  APPROXIMATION 

Introducing  the  g2  rules  the  constitutive  relations  for  the  fibers  and  matrices  become 


ME 

‘J 


-  Sfju  +  $!jET 
=  T 


(26) 

(27) 


The  effective  thermal  composite  properties  are  then  computed  as 

so  that 

PLASTIC  RESPONSE,  TRANSVERSELY  ISOTROPIC  APPROXIMATION 
Applying  the  g,  rules  to  equation  (13)  leads  to 

■MP  _  FMP±M.fMP  + 

zij  ~  Sijki  cki  +  Pi,  T 

so  that 

where 

fMT  _  aME  ,  fmp 
3 ijkl  -  ^ ijkl  +  ^ijkl 

and 

&r=&r+ftf 

The  effective  elastic-plastic  composite  properties  can  now  be  computed  as 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
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from  which  the  composite  elastic-plastic  constitutive  relation  becomes 


£ij  -  $ijki  Gy  +  Pj/  T 

(36) 

LAMINATE  MODEL 

The  constitutive  relations  (29)  and  (36)  can  now  be  used  to  build  an  incremental  laminate  model. 
In  order  to  do  so,  consider  a  laminate  of  thickness  2t  with  the  midthickness  reference  surface  located 
in  the  z=0  x-y  plane  of  the  right-handed  coordinate  system  x-y-z  (z  is  perpendicular  to  the  laminate 
layering).  Invoking  the  Love-Kirchoff  hypothesis  leads  to 

ux ( x ,y,z)  =  ux(x ,y)-z  ii°A (x ,y ) 

(37) 

Uy (x ,y ,z )  =  u°(x,y)-z  u°y (x  y ) 

(38) 

u2(x,y,z )  =  u°{x  j ) 

(39) 

where  the  functions  u°,  u7°,  are  the  displacements  of  the  reference  surface, 

corresponding  to  these  displacements  are 

The  strain  rates 

^xx  ~  ^xx  +  2  Kxx 

(40) 

Eyy  =  iyy+Z  Kyy 

(41) 

“  "Yxy  "*■ 2  Kx y 

(42) 

where  the  midsurface  strain  rates  are 

i°  —  u° 

cw  UXJC 

(43) 

i°  =  VL° 

yy  “y.y 

(44) 

+ 

v* 

ii 

(45) 
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and  the  rates  of  curvature  are 


**  = 


Kyy 


-z  u 


2  jcc. 


-Z  U, 


z.yy 


-  Z 


(46) 

(47) 

(48) 


Condensing  equations  (40),  (41),  and  (42)  wfth  matrix  notation  leads  to 


[£,]  =  [£*]  +  *[*] 


(49) 


Equations  (28)  and  (35)  can  also  be  expressed  in  terms  of  engineering  strains  and  matrix  notation  as 


[ei]  =  [S][a!]  +  [p]r 


(50) 


where,  depending  upon  the  state  of  the  layer,  [$)  and  [p]  are  either  elastic  or  elastic-plastic  matrices. 
Inverting  this  equation  leads  to 


[ail  =  iQ]Hi\+ [nr 


(51) 


where  [0]  is  the  layer  stiffness  matrix  and 


[H  =  -[£HP] 


(52) 


The  stress  and  strain  rates  expressed  in  the  local  layer  coordinate  system  can  be  transformed  to 
the  global  plate  coordinate  system  through 


[aj  =  [0o][a1] 

[£,]  =  [ec][e1] 


(53) 

(54) 
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where  the  transformation  matrices  are 


and 


ie0] 


cos20 
sin20 
sin0  cos0 


sin20  -2  sin0  cos0 

cos2©  2  sin0  cos0 
-sin0cos0  cos20-sin20 


[Be]  = 


cos20 

sin20 

2  sin0  cos0 


sin20  -sin0  cos0 
cos20  sin0  cos0 

-2sin0cos0  cos2©  -  sin20 


(55) 


(56) 


with  e  being  the  angular  orientation  between  the  local  1  and  the  global  x  coordinate  directions. 
Substituting  equations  (53)  and  (54)  into  equation  (51)  leads  to 

[<jJ  =  [G][ex]  +  [f]7  (57) 

where  [O]  and  [  r]  are  the  global  stiffness  and  thermal  stress  matrices  computed  from 

[e]  =  [e0][G][©£]"1  (58) 


and 


[n  =  [e0Hn  (59) 


From  the  laminate  edge  stress  rates  the  following  stress  resultant  rates  are  defined 

1 

(Nxx,Nyy,Nxy)=  l(pxx,cyy,6xy)dz  (60) 

-t 

t 

(Mxx,Myy,Mxy)=  jiOn'ayy'dxy)!  dz  (61) 

—t 

Switching  to  matrix  notation  and  introducing  equation  (57)  leads  to 


t  t  t 

[N]  =  j[Q]mdz  +  l[Q][K]zdz  +  j[T]Tdz  (62) 

-t  -i  -t 

i  i  t 

[M]  =  j[Q  ]  [e]  z  dz  +  j[Q ]  [K]  z2dz  +  J[f]  T  z  dz  (63) 

-i  -i  -i 
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which  can  be  combined  and  rewritten  in  the  familiar  laminated  plate  matrix  notation  as 


N 

A 

B 

i 

nt 

M 

__ 

B 

D 

k 

+ 

Mt 

(64) 


which  is  the  structural  stiffness  equation.  When  the  loading  rates  are  specified  this  equation  can  be 
inverted  to  yield  the  global  displacement  rates. 


PHASE  AVERAGE  MODEL 

After  equation  (64)  is  solved  the  laminate  equations  and  the  constitutive  relations  can  be  used  to 
determine  the  composite  stress  rates  for  each  layer.  The  constituent  phase  average  stresses  are  then 
computed  through  the  Phase  Average  model.  In  this  model  the  phase  average  stress  and  strain  rates 
are  related  through 


°ij  =  VF  + 

(65) 

*ij  -  VF  *ij  +  VM  *ij 

(66) 

where 


VM  =  1-VF  (67) 


Equations  (65)  and  (66)  are  a  consequence  of  using  the  CCA  as  the  micromechanical  model. 

For  layers  that  are  responding  elastically  the  phase  average  equations  can  be  combined  with  the 
constitutive  relations  to  yield 


For  layers  that  are  responding  plastically  equation  (68)  becomes 


These  equations  yield  the  matrix  phase  average  stress  rates.  The  fiber  phase  average  stress  rates  can 
then  be  determined  from  equation  (65). 
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APPLICATIONS 


DESCRIPTION 

The  engine  walls  of  the  NASP  vehicle  will  be  subjected  to  extremely  high  heat  fluxes.  To  devise  a 
structure  capable  of  withstanding  the  projected  thermal  loadings  It  has  been  proposed  that  the  walls  be 
designed  as  heat  exchangers  using  a  high  thermal  conductivity  material  that  possesses  high  stiffness 
and  strength.  Copper,  which  has  a  high  innate  conductivity,  has  its  stiffness,  strength,  and  conductivity 
enhanced  with  the  addition  of  advanced  Pl30x  graphite  fibers.  This  copper-based  MMC  will  be 
analyzed  here  assuming  a  fiber  volume  fraction  of  50%.  The  estimated  material  properties  for  the  Tiber 
and  matrix  are  listed  in  Tables  1  and  2. 

The  engine  walls  will  be  prechilled  prior  to  flight  from  an  assumed  ambient  atmospheric  temperature 
of  20  deg  C  to  the  fuel  temperature  of  -170  deg  C.  During  flight  the  temperature  of  the  engine  walls 
are  anticipated  to  reach  650  deg  C.  After  the  flight  the  walls  will  return  to  an  ambient  condition.  In  this 
work  the  described  thermal  cycle  is  analytically  applied  to  both  unidirectional  and  cross-ply 
Pl30x/copper  laminates. 

NUMERICAL  IMPLEMENTATION 

The  theoretical  model  was  developed  in  terms  of  rates.  However,  the  application  of  the  model  is 
necessarily  incremental  with  piecewise  linearity  assumed  for  each  step  of  the  analysis.  This  is  an 
approximation  since  the  system  characteristics  are  temperature  and  stress  dependent  and  can  vary 
continuously  over  the  course  of  a  load  step.  To  correct  for  this  and  for  the  approximations  of  the 
transversely  isotropic  assumptions,  the  numerical  model  iterates  to  predict  the  onset  of  yield  and  to 
determine  the  extent  of  plastic  flow.  An  evaluation  of  the  effectiveness  of  these  numerical  procedures 
can  be  found  in  reference  13. 

An  additional  control  over  the  accuracy  of  the  analysis  is  the  size  of  the  load  increments.  Table  3 
shows  the  results  of  a  numerical  study  in  which  successively  smaller  temperature  increments  where 
used  in  the  analysis  of  a  unidirectional  Pl30x/copper  plate  subjected  to  one  thermal  cycle.  A 
continuous  plot  of  the  response  of  the  plate  for  half-degree  temperature  increments  is  shown  in  Figure 
1.  (An  explanation  of  the  predicted  response  can  be  found  in  the  following  subsection, 
“P130x/Copper*.)  For  temperature-dependent  hardening  the  field  variables  converge  with  decreasing 
temperature  increments  while,  as  a  consequence  of  the  micromechanical  model,  the  composite  and 
matrix  axial  strains  align. 

Table  3  and  Figure  1  also  show  the  results  of  a  temperature-independent  hardening  model.  Note 
that  this  assumption  introduces  a  slight  deviation  from  the  more  general  model  at  the  high  end  of  the 
temperature  cycle  but  achieves  convergence  of  the  field  variables  at  a  larger  temperature  increment. 

P130X/COPPER 

The  first  problem  to  be  considered  is  a  unidirectional  plate  subjected  to  two  successive  thermal  load 
cycles.  Figure  2  is  a  plot  of  the  composite  axial  strain  for  this  structure  in  which  the  response  is  seen 
to  be  extremely  nonlinear  and  hysteretic  but  with  no  ratcheting.  The  composite  initially  responds  during 
precooling  by  contracting  elastically.  However,  the  matrix  material  yields  at  -23.2  deg  C  after  which 
the  fiber  properties  dominate.  Hence,  with  further  cooling  the  composite  begins  to  grow  axially.  On 
heating,  the  matrix  response  is  again  elastic  so  that  the  composite  continues  to  grow  untB  yield  occurs 
at  -61 .5  deg  C.  From  here  on  the  composite  contracts  untl  the  changing  axial  and  transverse  material 
properties  interact  to  produce  thermal  growth.  At  the  high  end  of  the  thermal  cycle  the  growth  peaks. 
With  cooling,  the  composite  contracts  elastically  and  then  plastically  (yielding  at  630  deg  C).  Further 
cooling  leads  to  axial  growth.  The  next  and  successive  thermal  cycles  lead  to  an  identical  response. 
Figure  3  plots  the  composite  transverse  strain  in  which  it  is  seen  that  the  nonlinearity  is  not  as  severe 
since  the  fiber  and  matrix  thermal  expansivities  are  of  the  same  sign.  A  slight  hysteresis  is  also  evident 
here.  Figures  4  and  5  plot  the  matrix  plastic  strains  which  occur  isovolumetrically.  Figures  6  and  7 
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Table  1 .  Fiber  Material  Properties. 


P130x  Graphite  Fiber 

Temperature 

a* 

a, 

deg  C 

\izper  deg  C 

pe  per  deg  C 

-200 

-1.62 

10.0 

-100 

-1.62 

10.0 

0 

-1.62 

10.0 

20 

-1.62 

10.0 

100 

-1.60 

10.8 

200 

-1.05 

10.8 

300 

-.625 

11.6 

400 

-.275 

12.4 

500 

.100 

12.6 

600 

.450 

12.6 

700 

.815 

_ _ 

12.6 

aa 

axial  secant  thermal  expansion  (reference  temperature  =  20  deg  C) 

at 

transverse  secant  thermal  expansion 

Ea  =  938  GPa, 

axial  extensional  modulus 

£f  =  19.3  GPa, 

transverse  extensional  modulus 

va  =  0.210, 

axial  Poisson's  ratio 

v(  =  0.892, 

transverse  Poisson's  ratio 

Ga  =  17.9  GPa, 

axial  shear  modulus 
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Table  2.  Matrix  Material  Properties. 


Copper,  OFHC  Grade 

Temperature 

£ 

V 

a 

£ 

deg  C 

GPa 

peper  deg  C 

MPa 

Gpa 

-200 

137 

.360 

14.4 

107. 

8.77 

-100 

126 

.355 

15.3 

86.4 

8.29 

0 

118 

.345 

16.3 

72.5 

8.05 

20 

117 

.343 

16.5 

70.0 

8.00 

100 

114 

.333 

17.0 

60.5 

7.58 

200 

110 

.322 

17.4 

52.0 

6.94 

300 

106 

.312 

17.7 

42.6 

6.12 

400 

102 

.307 

18.1 

31.5 

5.14 

500 

98 

.299 

18.5 

21.4 

4.56 

600 

93 

.288 

19.0 

15.1 

4.00 

700 

88 

.275 

19.5 

11.7 

3.44 

E  extensions  modulus 

v,  Poisson's  ratio 

a  secant  thermal  expansion  (reference  temperature  «  20  deg  C) 

oy  yield  stress 

tangent  modulus 


14 


NADC-90004-60 


Table  3.  Numerical  Study  of  a  Pl30x/Copper  Unidirectional  Plate. 


Temperature  Dependent  Hardening  (1)  (2) 

Step  Size 

CPU  Time 

-170  deg  C 

650  deg  C 

Cool  Down 

deg  C 

Sec 

mm 

MhM 

EH 

mm 

WEM 

eh 

msm 

tii 

Bl 

5.0 

51 

1.358 

1.321 

18.83 

6.022 

5.314 

-97.60 

-.8899 

-2.166 

-7.309 

2.0 

119 

1.359 

1.344 

18.82 

6.089 

5.806 

■9122 

-.7613 

-1.272 

1.0 

217 

1.359 

1.352 

18.82 

6.111 

5.970 

-97.10 

-.7185 

-.9737 

-6.378 

0.5 

425 

1.359 

1.355 

18.82 

6.123 

6.052 

-97.03 

-.6970 

-.8247 

-6.261 

Temperature  Independent  Hardening  (1)  (2) 

Step  Size 

CPU  Time 

-170  deg  C 

650  deg  C 

Cool  Down 

deg  C 

Sec 

worn 

^11 

EH 

mm 

eff 

mm 

■EH 

7* 

en 

mm 

5.0 

49 

1.369 

1.369 

18.% 

6.444 

6.444 

-92.07 

-.8999 

-.8999 

-6.143 

2.0 

I- 

112 

1.370 

1.370 

18.93 

6.512 

6.512 

-92.09 

-.7653 

-.7653 

1.0 

213 

1.370 

1.370 

18.92 

6.535 

6.535 

-92.10 

-.7205 

-.7205 

-6.144 

0.5 

418 

1.371 

1.371 

18.92 

6.547 

6.547 

-92.10 

-.6980 

-.7018 

-6.148 

(1)  The  exponent  of  the  strains  is  -4. 

(2)  The  computations  were  performed  on  a  Sun  Microsystem  3/280. 


15 


P130x/Copper,  Unidirectional 


NADC-90004-60 


urejjs  ejjsodujoo  |B|XV 


16 


Figure  1.  Composite  Axial  Strain  vs.  Temperature  for  Both  Temperature- 
Dependent  and  Temperature-Independent  Hardening. 
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show  the  matrix  stresses  caused  by  the  mismatch  in  constituent  straining.  Note  that  in  this  problem  the 
thickness  matrix  stress  equals  the  transverse  stress,  while  the  shear  stresses  are  zero. 

Figure  8  plots  the  load  function  over  the  thermal  cycle.  This  figure  is  useful  in  interpreting  the 
different  regions  of  thermoplastic  response.  During  the  initial  cooling  the  matrix  is  elastically  restrained 
by  the  fibers  so  that  the  load  function  quickly  reaches  the  yield  line  along  which  it  travels  plastically. 
(The  yield  line  is  a  plot  of  the  strength  term  of  equation  (9).)  Heating  from  the  low  temperature  causes 
the  load  function  to  unload  and  reload  elastically  and  then  travel  along  the  yield  line  plastically.  In 
stress  space  this  is  equivalent  to  unloading  from  the  yield  surface,  transversing  the  elastic  region,  and 
then  striking  the  yield  surface  at  another  point.  The  same  behavior  occurs  when  the  plate  is  cooled  from 
the  high  temperature  point.  Subsequent  cycles  lead  to  identical  results.  Note  that  the  effect  of  the 
temperature  dependency  is  reflected  in  the  length  of  the  elastic  path  at  the  two  temperature  extremes. 

Because  of  the  method  of  manufacture  MMC  typically  contains  residual  processing  stresses. 
Although  intuition  would  lead  one  to  suspect  that  these  stresses  are  tensile,  there  is  evidence  showing 
that  surface  friction  and  high  processing  pressures  lead  to  residual  stresses  that  are  actually 
compressive.  The  second  problem  to  be  considered  analyzes  the  effect  of  residual  processing  stresses 
on  the  thermal  response  by  assuming  a  residual  stress  state  of 

(off,  Oj2,  ag.  og)  =  (-80,-40,-40, 0, 0, 0)  MPa 

in  a  unidirectional  plate.  Figure  9  is  a  plot  of  the  load  function  of  this  plate  when  it  is  subjected  to  two 
thermal  cycles.  The  path  starts  at  the  residual  stress  state,  unloads,  reloads  and  then  reaches  the  yield 
line  at  -60.8  deg  C.  The  subsequent  load  path  and  the  location  of  the  yield  points  are  virtually  identical 
to  the  case  with  no  residual  processing  stresses  (Figure  8).  Because  of  this  the  shape  of  the  curves  of 
the  field  variable  in  the  residually  stressed  plate  (Figures  10  through  15)  are  very  similar  to  the  plots  of 
the  initially  unstressed  plate.  The  difference  in  response  is  seen  to  be  a  vertical  shift  of  the  curves  and 
hence  the  peak  magnitudes,  if  a  residual  tensile  stress  state  had  been  examined,  then  similar  shifting 
would  occur  except  in  an  opposite  direction. 

The  next  problem  to  be  considered  is  a  0/90/90/0  cross-ply  laminate  (no  residual  stresses). 
Figure  16  shows  that  for  such  a  structure  the  composite  axial  strains  will  ratchet  under  cyclic  thermal 
loads.  Because  the  construction  is  cross-plied  and  symmetrical,  the  transverse  strains  are  identical  to 
the  axial  strains.  The  matrix  plastic  strains  (Figures  17  through  19)  show  a  pronounced  ratcheting. 
Since  the  construction  constrains  the  plys  from  deforming  individually,  the  thermal  cycle  will  result  in 
non-zero  composite  stresses.  Figures  20  and  21  plot  these  composite  stresses  for  the  upper  0-deg  ply. 
The  matrix  stresses  are  therefore  a  superposition  of  the  effects  of  global  composite  stressing  and  of 
local  differential  thermal  growth.  Figures  22  through  24  show  that  these  stresses  ratchet  substantially 
during  the  thermal  cycle.  Without  the  processes  of  stress  relaxation,  unacceptably  high  stress  levels 
would  soon  be  reached. 
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Figure  2.  Composite  Axial  Strain  vs.  Temperature. 


leuojpejipjun  ‘jaddoo/xoe  l-d 


NADC-90004-60 


UIOJ^S  ©SJ&ASUDJJ.  0}ISOdLUO£) 


19 


Figure  3.  Composite  Transverse  Strain  vs.  Temperature. 
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•e  5.  Matrix  Transverse  Plastic  Strain  vs.  Temperature. 
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Figure  8.  Load  Function  vs.  Temperature. 
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Figure  9.  Load  Function  vs.  Temperature. 
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Figure  10.  Composite  Axial  Strain  vs.  Temperature. 
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Figure  11.  Composite  Transverse  Strain  vs.  Temperature. 
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Plastic  Strain  vs.  Temperature. 
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Figure  13.  Matrix  Transverse  Plastic  Strain  vs.  Temperature. 
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Figure  14.  Matrix  Axial  Stress  vs.  Temperature. 
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Figure  15.  Matrix  Transverse  Stress  vs.  Temperature. 


0/06/06/0  ‘Jaddoo/xoetd 


NADC-90004-60 


uioj^s  (Dixy  ejisoduuoo 


32 


Figure  16.  Composite  Axial  Strain  vs.  Temperature. 
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Figure  18.  Matrix  Transverse  Plastic  Strain  vs.  Temperature. 
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Thickness  Plastic  Strain  vs.  Temperature. 
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Figure  21.  Composite  Transverse  Stress  vs.  Temperature. 
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Figure  23.  Matrix  Transverse  Stress  vs.  Temperature. 
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CONCLUSIONS 

In  this  work  a  thermoplasticity  model  was  formulated  for  the  purpose  of  examining  the  mechanical 
and  thermal  behavior  of  MMC.  The  theoretical  model  was  applied  to  the  analyses  of  materials  that  are 
currently  being  proposed  for  use  on  advanced  aerospace  vehicles.  The  applications  demonstrate  the 
usefulness  of  the  theory  in  predicting  and  interpreting  basic  material  response.  Also,  the  analysis  data 
can  be  quickly  changed  to  parametrically  and  inexpensively  assess  other  designs. 

For  the  Pl30x/copper  material  the  analysis  predicts  that  thermal  cycling  will  lead  to  substantial 
plastic  straining  and  hysteresis.  Also,  depending  upon  the  construction,  the  hysteresis  loops  can  ratchet 
under  successive  thermal  loadings.  The  effect  of  the  residual  processing  stresses  was  seen  to  be  an 
upward  or  downward  shifting  of  the  response  curves. 
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